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COMPLEX PERMITTIVITY MEASUREMENTS DURING HIGH TEMPERATURE RECYCLING OF 


SPACE SHUTTLE ANTENNA WINDOW AND DIELECTRIC HEAT SHIELD MATERIALS 

By Harold L. Bassett and Steve H. Bomar, Jr. 

Engineering Experiment Station 
Georgia Institute of Technology 

SUMMARY 

Thermal and electrical measurements have been performed on six candidate 
space shuttle antenna window materials: Lockheed LI-1500, McDonnell Mullite 

HCF, hot pressed boron nitride, Philco-Ford AS-3DX, General Electric Markite 
3DQ, and slip-cast fused silica. The material samples were cycled ten times 
through a simulated reentry surface temperature profile. 

During temperature cycling, surface and internal temperatures were 
measured and the complex microwave transmission coefficient was obtained 
using a free-space phase shift bridge technique. Using these measured data, 
the dielectric constant and loss tangent of each material as functions of 
temperature were calculated to determine whether these properties deterio- 
rated with repeated temperature cycling. 



INTRODUCTION 


Scope of the Program 

The scope of this research program was to perform measurements to 
determine the values of dielectric constant and loss tangent of several 
candidate Space Shuttle antenna window and heat shield materials as detailed 
functions of temperature during exposure to high temperature recycling tests 

Modifications to the high temperature complex permittivity measurement 
system developed by the Engineering Experiment Station of the Georgia 
Institute of Technology under two Air Force Contracts (Ref. 1,2) were made 
so that the complex permittivity of candidate heat shield and antenna window 
materials could be measured during exposure to Space Shuttle reentry surface 
temperature profile recycling tests. 

Statement of the Problem 

The change in electrical properties of electromagnetic window materials 
as a function of temperature recycling is not known. Since the basic struc- 
ture of the NASA Space Shuttle will be used on ten or more flights, it was 
believed necessary to determine the change in electrical properties of 
candidate electromagnetic window materials as a function of reentry flight 
time and as a function of temperature recycling. Each candidate material 
was subjected to a temperature versus time profile designed to simulate 
the anticipated reentry flight profile. The materials were heated on one 
side by gas torches and the rate of heating was controlled during the 
measurement runs to obtain the required surface temperatures on the samples. 
Each sample was subjected to ten temperature cycles. Two samples of each 
material type were run. A free -space microwave phase shift bridge was used 
to measure the change in transmission properties of the materials as a func- 
tion of time. Dielectric constant and loss tangent data were obtained from 
the transmission and temperature data. 

Background 

There has been considerable research in the determination of the elec- 
trical properties of materials as functions of temperature. Of particular 
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interest are the works of Westphal (Ref* 3,4), Bowie (Ref. 5), and Gilreath 
(Ref. 6). These researchers have used the short circuited waveguide tech- 
nique (Ref. 7,8) to determine the complex permittivity of materials as a 
function of temperature. In the short circuited waveguide technique the 
material of interest is inserted into a metallic sample holder which is 
enclosed within a furnace. Microwave transmission measurements are then 
made to determine the complex permittivity. Excellent results are obtained 
from this technique to temperatures approaching 1923° K (3000° F) . The 
material sample is heated in a near uniform manner so that the whole sample 
is tested at the particular temperature of interest. 

For the recycling measurements it was felt necessary to subject the 
samples to the anticipated flight thermal environment. This necessarily 
meant that one surface of the samples would be exposed to the heat source. 
Then, to obtain useful data, a free-space microwave measurement system was 
employed in obtaining the microwave transmission properties of the material 
sample as a function of run time. 
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MEASUREMENT SYSTEM 


The focused beam system allows free-space phase and attenuation measure- 
ments to be made on electromagnetic window materials at X-band frequencies 
to very high temperatures. In fact, measurements have been made to temper- 
atures above 2480° K (4000° F) . For this particular program the peak tem- 
perature was about 1538° K (2300° F) . The complex permittivity versus tem- 
perature function is determined from the transmission properties of the 
material. Past research programs at Georgia Tech have utilized this system 
at both X-band and Riband frequencies (Ref. 1,2). 

Description of Free-Space Microwave Sysgem 

The complete microwave system is diagrammed in Figure 1. The basic 
system employs the free space phase-shift bridge and transmission technique 
described by Redheffer (Ref. 9) to determine the complex permittivities of 
materials. The focused beam system forms the sample arm of the microwave 
bridge. The other arm of the bridge is the reference path for the microwave, 
signal. The material specimen is placed in the focal region of the prolate 
spheroidal focusing reflector normal to the beam axis in the sample arm, and 
the resulting change in electrical -path length is measured by an audio phase 
detector network. Simultaneously, the attenuation due to the sample is 
measured. The complex permittivity is computed from these transmission data. 
Electrical measurements can be performed either in a static or a dynamic 
thermal environment since the phase and amplitude data are measured and 
recorded on a continuous basis without having to adjust a mechanical phase 
shifter. 

The technique used to obtain the phase and amplitude data is as follows: 
A phase-locked CW carrier frequency, f Q , is used both as the frequency of 
transmission for the focused beam and for the generation of a reference 
signal, f + 1000 Hz. As depicted in Figure 1, the reference signal is 
mixed with the carrier signal at three single-ended crystal mixers to obtain 
incident, reflected and transmitted signal measurements. The 1000 Hz mixer 
outputs are amplified and recorded for amplitude data. Simultaneously, the 
phase characteristics 0^ and of the transmitted and incident signals 
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respectively are measured by two phase detectors. The 1000 Hz input reference 
signal to the phase meters is the same signal that is used to modulate the 
single-siaeband generator network. These phase and amplitude data are 
recorded on a multichannel recorder. 

The single-sideband generator utilizes a 3-dB sidewall hybrid coupler 
terminated with PIN diodes modulated by a 1000 Hz square-wave signal to 
generate the + 1000 Hz reference signal. The slide screw tuners, as 
shown in Figure 1, are adjusted to obtain the proper amplitude and phase 
relationships of the signal in the two hybrid arms. 

The microwave system, as described, can be used to measure the trans- 
mission properties of a material at (1) ambient temperature, (2) uniform 
elevated temperature, and (3) in a transient state (temperature gradient). 

The latter measurements are made as a function of time and correlated with 
temperature data to calculate the complex permittivity of the material as a 
function of temperature. The measurements are made at a frequency of 9.313 
GHz. 

The focused beam is formed by a 117 cm (46-inch) diameter prolate 
spheroidal reflector which is illuminated by a feed horn placed at the inner 
focal point of the reflector. The material sample is positioned at the outer 
focal point of the reflector which is located 109 cm (43 inches) from the 
reflector apex. The incident beam is focused to a 3.18 cm (1.25-inch) 
diameter circle (3 dB points) at the outer focal point. An identical 
reflector, whose outer focal point is coincident with that of the trans- 
mitting reflector, acts as the receiving aperture in the free-space bridge. 

The sample is mounted on a 51 cm (20.5 inch) diameter carrier disk which is 
enclosed in a furnace that has open windows on each side through which the 
microwave beam passes. The disk-shaped sample carrier is mounted so that 
it is revolved about its center, and the microwave beam is Incident on a 
circular area whose center is displaced from the disk center by 17.1 cm 
(6-3/4 inches). 


Description of Sample Rotating 
and Heating Apparatus 

Samples were mounted in a 51 cm (20.5 inch) diameter carrier disc 
which rotated inside a furnace enclosure. The rotating carrier and mounted 
samples were heated by natural gas-air flames using techniques very similar 
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to those developed on two Air Force programs which preceded this contract 
(Ref. 1 and 2). The only substantial modifications of earlier experimental 
techniques involved heating the samples for long periods of time through a 
predetermined temperature versus time profile. Formerly, emphasis had been 
placed on achieving the highest possible temperatures without regard to the 
time the run was in progress. 

Design and Construction of Apparatus . The microwave waveguides and 
reflectors were in a fixed position and could not easily be moved. Thus, 
the sample rotating and heating device was designed to roll on tracks; 
movement in two directions was provided. First the entire furnace could 
be translated parallel to the microwave beam to allow measurement of the 
average value of transmission coefficient. During measurement runs, the 
furnace was kept in continuous motion in this direction through a distance 
of i 1 centimeter from the center position. The second direction of motion 
was perpendicular to the microwave beam, for installation of samples, cali- 
bration and servicing of the equipment. 

Figures 2 and 3 show two views of the sample heating furnace and other 
parts of the system. The basic frame was constructed in a configuration 
somewhat like a sawhorse on wheels. The mandrel on which the samples turned 
was located at the end of a horizontal arm which was pivoted at the top of 
the sawhorse frame to allow vertical adjustment of the sample position. The 
front half of the furnace box pivoted upward as shown in Figure 3 to permit 
sample installation and other servicing operations. An electric motor, 
mounted near the bottom of the frame, drove the rotating mandrel through a 
speed-reducing gear box and chain drive. A carrier disc with several mounted 
samples is installed in the furnace and can be seen in the figures. 

On the raised section of the furnace, the cone which defines the micro- 
wave beam and the gas torch plumbing are visible. Furnace walls were lined 
with 2.5 cm (1 inch) thick fused silica foam refractory blocks, and the 
cone was cast from fused silica and lined on the heated side with silica 
foam. The torches were made up at Georgia Tech using commercially available 
cast iron burner heads and rigid copper tubing. The burner heads were water 
cooled; natural gas and air were supplied from manifolds near their respec- 
tive control valves, and the gas and air were mixed in the large tubes 
leading to the burner heads. Supplementary oxygen was injected into the 
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combustion air at certain times during the runs to increase combustion 
temperatures. The six gas torches were located at radii on the carrier 
which were selected to give uniform heating across the faces of the samples. 

The temperature was controlled by manual adjustment of gas, air, and 
supplementary oxygen flow rates. Time and funds did not permit installation 
of an automatic temperature control system. The "control temperature" was 
measured by a thermocouple located on the heated carrier surface. The out- 
put of this thermocouple was plotted on an x-y recorder; prior to the run 
a curve representing the desired front surface temperature versus time was 
plotted on the recorder chart, so that the furnace operator had a visible 
record of these temperatures to aid in controlling the furnace. 

Samples were centered on a circle at a radius of 17 cm (6.75 in.) from 
the center of the carrier disc. The furnace was positioned so that the 
microwave' beam passed through the cone, the samples, and the beam hole in 
the rear furnace wall. When the system was in operation no part of the 
furnace intruded into the microwave beam. The beam occupies roughly a 
cone shaped volume in the region near the sample, so that the flame control 
baffle must be cone shaped. 

It was necessary to exercise care to prevent flame attenuation of the 
microwave beam. In this program, the gas, air, and oxygen flow rates to 
the torches were adjusted throughout the runs. The ratio of fuel to oxi- 
dizing gases was the main factor controlling flame interference; too much 
fuel gave visible flame in the microwave beam area and significant attenuation 
could be observed on the microwave recording instruments when this condition 
existed. Thus, conscious efforts were made to assure that the torches 
operated with an oxidizing flame. A hood was installed above the sample 
rotating and heating apparatus for removal of combustion gases from the 
room, but in contrast to earlier programs this exhaust system played no part 
in flame control. 

The carrier disc was supported on the mandrel assembly shown in Figure 
4. The large bearing blocks contained channels for cooling water, and water 
was continuously circulated through the mandrel assembly. The rotating 
shaft was provided with 18 small conduits through which thermocouple wires 
were passed from the rear of the carrier disc to brass slip rings outside the 
furnace. The slip rings rotated in small troughs filled with mercury to 
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conduct thermocouple signals from the mandrel to the recorders. 

Two lock nuts were used to retain the carrier disc in place on the end 
of the mandrel and the nuts were covered by a block of silica foam to protect 
the threads from damage by flame; also the threaded stud was hollow to permit 
water cooling. A microswitch and cam were installed on the shaft adjacent 
to the slip rings to provide an event mark in each revolution; these event 
marks were placed on all recordings. 

Techniques of Surface Temperature Measurement . Optical measurement 
surface temperatures had been accomplished on an earlier program by the use 
of fast-response infrared pyrometers. The instruments were Barnes Engineer- 
ing Company Model IT-4 pyrometers with several custom modifications: 


(1) Narrow-band optical filters and un immersed thermistor bolometers 
(detectors) ; 

(2) Field of view subtending a 2 cm (0.75-in. ) diameter spot at a 
range of 1.8 m (6 ft.); 

(3) Output range from ambient temperature to 2750° K (4500° F) without 
range switches; and 

(4) Electrical output 0 to 5 volts over the temperature range of the 
instrument with a time constant of 15 milliseconds. 


The pyrometers were intended to have sufficiently fast response that 
the temperatures of individual samples could be measured while the samples 
were passing by the microwave window. Also, they operated at spectral 
wavelengths where silica-based materials have high emissivity; passbands of 
4.5 pm, 5.1 p m, and 7.9 pm could be obtained with the available filters. 

The sensing head of one of these pyrometers is visible in Figure 2, 
mounted on a stand in front of the right microwave reflector. The pyrometers 
were sensitive to high ambient temperatures, so that a cooling shield was 
provided and the shield was cooled by air from a small air-conditioning unit. 
The field of view was much larger than specified, so that the head required 
mounting at a short stand-off distance in front of the microwave reflector. 

The pyrometer output voltage is directly proportional to target radiance. 
For a narrow band of wavelengths, target radiance is related to temperature 
by Planck's law (Ref. 10): 


-5 


b,A 


= 2 C 


1 C./XT 
e 2 


( 1 ) 
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where 


N, = radiant intensity at X 

o > A 

-13 2 

= first constant of Planck's Law = 5.956 x 10 watt cm 

4 « 

= second constant of Planck's Law = 1.4387 x 10 IJm K 


X = wavelength in pm 
and T = black body temperature in °K. 

For a band of wavelengths 




N b = 


N b,l dX 


( 2 ) 


where 

= radiance = radiant energy per unit of time, 
unit of projected area and unit of solid angle. 

This function is not easily integrated between limits X^ and . However, 
the bandpass interval is very narrow in this case so that we may use the 
approximation : 


N b = N b,X AX = 


2 C x A X 

X 5 (e C 2 /XT -l) 


(3) 


Barnes Engineering Company furnished a table of radiance vs temperature 
for each filter and also the cutoff and peak values of wavelength for each 
filter. Thus sufficient data were available to permit calculation of 
surface temperatures from the recorded pyrometer output if a calibration 
could be established. Since the instrument output is a linear function of 
target radiance, and radiance defines the black-body temperature through 
Planck's Law, the calibration step must relate output voltage to target 
radiance. Calibration was performed by measuring the instrument outputs 
when the pyrometer was pointed at a slip-cast fused silica plate whose 
temperature was known. The plate was placed behind the door of an electri- 
cally heated furnace. The furnace temperature was allowed to equilibrate and 
the temperature of the target was measured by a thermocouple embedded under 
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its surface. The pyrometer was pointed at the furnace door and the recorder 
started. The furnace door was opened quickly for a few seconds to obtain a 
recording of the pyrometer's response to the hot target. The target tem- 
perature began to drop immediately upon opening of the furnace door; thus 
the pyrometer output at the instant of opening was used for calibration. 
Calibration data were obtained at several temperatures over the range of 
interest on this program, and a plot of output voltage vs target radiance 
was a straight line as expected. 

The optical pyrometer recordings were converted to front surface 
temperatures by a program written for an electronic calculator. The inte- 
grated Planck equation (equation 3) can be rearranged to express the tem- 
perature as : 

m C 2 (4) 

2C AX 

X In [1 + — ] 

x\ 

The program was written to input the millimeters of pen deflection read 
from the recording paper, calculate from the calibration line, calculate 
the temperature in °K, and print the temperature in °F. A separate program 
was prepared for each pyrometer, incorporating the appropriate calibration 
and wavelength constants. 

The front surface temperatures reached on this program were low enough 
to permit the use of front surface thermocouples. Thus, platinumrplatinum- 
13% rhodium couples were installed on the front surface of each thermal 
sample to check the optical data. It is recognized that metallic thermo- 
couples do not have the same surface characteristics as the ceramic sample 
materials, and that their temperature responses might not therefore be 
exactly representative of the samples. However, generally satisfactory 
agreement between optical and thermocouple temperatures was obtained. 
Thermocouples have the advantage of better resolution at low temperatures 
than optical pyrometers. 


14 



preparation of samples 


Six sample materials, candidates for antenna windows in the Space Shuttle 
thermal protection system, were investigated on this program. Two sets of 
samples for each material were processed through the test sequence. A set 
of samples consisted of a nominal 15.2 cm (6-in.) diameter microwave sample 
and a nominal 3.8 cm (1.5-in.) diameter thermal sample; all thermocouples 
were mounted in thermal samples because any metal in the microwave samples 
would interfere with radar measurements. 

The preferred sample thickness for these measurements is a multiple of 
one-half wavelength at room temperature. Adherence to this starting thick- 
ness reduces reflections in the microwave system. However, low density ma- 
terials were run at the thickness as received because trimming would have 
reduced the quantity of material per unit area and thereby reduced the pre- 
cision of microwave phase shift and attenuation measurements. 

A list of the sample materials and their manufacturers is given in Table I. 
All samples were supplied to Georgia Tech by NAS A- Lang ley except slip-cast 
fused silica. 


Mounting of Samples 

The microwave and thermal samples were mounted in fused silica carrier 
disks as shown in Figure 5. The carriers were made by pouring a wet casting 
mixture of fused silica grain and slip into a gypsum plaster mold; the mixture 
sets as water is drawn into the plaster mold. Holes for the mandrel and 
samples were formed by locating dense fiberboard plugs of appropriate sizes 
on the mold before pouring in the casting mixture. After the carrier had 
set up, the hole plugs were removed and the carrier was dried for several 
days, then fired at about 1475° K (2200° F) to develop strength. 

Chrome 1 - alume 1 embedded and backside thermocouples and platinum-rhodium 
front side thermocouples were installed in the thermal samples before samples 
were mounted in the carriers. Thermocouple leads were made long enough to 
reach the slip rings after the carrier was mounted on the mandrel. Two 
embedded thermocouples were positioned along a diameter of the sample, in a 
1.2 mm (0.030 in.) wide saw cut at the desired depth. Saw cuts were made 
from the rear (unheated) side of the sample, and successive cuts were placed 
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TABLE I 


SAMPLE MATERLALS AND MANUFACTURERS 


Description 

Nomenclature 

Manufacturer 

Fused quartz reinforced 

AS-3DX 

Phi 1 c o - Ford Corpora t ion 

silica compo site 



Slip- cast fused silica 

SCFS 

Georgia Tech 

Fiber bundle reinforced 

Markite 3-DQ 

General Electric Co, 

silica composite 



Silica fiber rigidized 
coating 

LI-1500 

Lockheed Missiles and 
Space Co, 

Rigidized mullite fiber 

Mullite HCF 

Me - Donne 1 1 - Doug la s 
Astronautics Co . 

Hot pressed boron nitride 

HD-0092 

Union Carbide Corporation 
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Figure 5. Front View of Carrier Disk and Samples. 



on different diameters so that the wires were parallel to the plane of the 
sample surface but not parallel in other planes. After the wire was positioned 
in the cut, the cut was filled with a cement made of Ludox AS* (colloidal 
silica) and powdered sample material. The saw cuts were carefully filled to 
exclude air bubbles, and the depth of each thermocouple was recorded. After 
experimental runs, the thermal samples were sectioned to verify thermocouple 
depths. 

The samples were inserted in their respective holes in the carrier disk 
and anchored by flowing silica cement into the space between the sample and 
carrier. The samples were positioned with their heated surfaces flush with 
the front surface of the carrier, as shown in Figure 5. A photograph of the 
back of this same carrier is shown in Figure 6. Thermocouple wires were covered 
with silica cement across the back of the carrier to prevent their being damaged 
by flames during the ten temperature cycles. 

Physical Properties of Sample Materials 

The physical data given in Table II are for the purpose of documenting 
the several sample materials. This list is not intended to show all available 
physical data, but rather to provide data for identifying and characterizing 
the samples. 

Mullite HCF and LI-1500 were run with coatings on the heated surface of 
the samples. These emittance control coatings were applied by the manufac- 
turers of the two materials, and their compositions and properties were not 
known to Georgia Tech investigators. In each case, the coatings were less 
than one millimeter thick. 

Carrier and Sample Assignments 

Four carriers were run to obtain the required measurements after assembly 
and check out of the system had been completed. Duplicate specimens of 
each material were carried ten times through the simulated reentry temperature 
versus time profile. Samples were assigned to carriers in a manner that pre- 
cluded running both samples of one material on the same carrier. The materials 
run on each carrier are shown in Table III. 


Trade Mark of E. I. duPont de Nemours and Co. 
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Figure 6. Rear View of Carrier Disk and Samples. 



TABLE II 


PHYSICAL PROPERTIES OF SAMPLE MATERIALS 


Slip-Cast Fused Silica (SCFS) 

Slip; Thermo -Mate rials Corporation, batch 041772-2, high purity 

Slip Median Particle Size: 6.9 pm 

Firing Conditions: 1478° K (2200° F) , 4 hours 

Bulk Density: 1.902 gm/cm^ 

Silica Content: 99.56 weight per cent 

Hot Pressed Boron Nitride (HD-0092) 

Boron content: 43.0 to 43.5 per cent 

Nitrogen content: 55.6 to 56.3 per cent 

Oxygen content: < 1.00 per cent 

Carbon content: <0.4 per cent 

Other metallic impurities: 0.05 per cent 

Bulk Density: 1.945 gm/cm^ 

Fused Quartz Reinforced Silica Composite (AS-3DX) 

3 

Bulk Density: 1.673 gm/cm 

Fiber Bundle Reinforced Silica Composite (Markite 3DQ) 

3 

Bulk Density: 1.913 gm/cm 

Quartz Fibers: 61.47 weight per cent 

Silica: 38.53 weight per cent 

Porosity: 11.81 volume per cent 

Silica Fiber Rigidized Coating (LI-1500) 

3 

Bulk Density: 0,237 gm/cm (uncoated) 

Rigidized Mullite Fiber (Mullite HCF) 

3 

Bulk Density: 0.279 gm/cm (uncoated) 

Porosity: 86 volume per cent 
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TABLE III 


CARRIER AND SAMPLE ASSIGNMENTS 


Carrier Number 

Samples 

i 

SCFS , Markite 3DQ, Boron Nitride HD-0092 

2 

AS-3DX, Mullite HCF, LI-1500 

3 

Boron Nitride HD-0092, AS-3DX, Mullite HCF 

4 

SCFS, Markite 3DQ, LI- 1500 
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EXPERIMENTAL MEASUREMENT PROCEDURE 


The carrier disk containing the mounted samples is installed on the 
sample rotating and heating apparatus, and thermocouple connections to the 
recorders are completed. The microwave equipment is allowed to warm up for 
at least one hour prior to calibration of the system. 

To begin calibration, the klystron signal frequency is adjusted to 
9.313 GHz and is phased- locked by a synchronizer which maintains a frequency 
of 9.313 GHz ± 0.0007 MHz. E-H tuners are used to reduce the Voltage Standing 
Wave Ratio to less than 1.02 at the output and input of the measurement arm, 
with no furnace or other obstruction in the beam. (It has previously been 
established that an empty furnace located in the measurement position has no 
influence on the microwave beam.) 

The free-space phase and amplitude calibrations are then made. The 
incident and transmitted signal calibrations are made by an adjustment of the 
precision variable attenuator which is shown in Figure 1. The phase calibra- 
tions of the transmitted signal are made by adjusting the phase shifter shown 
in Figure 1. During this operation, signals having known values are placed 
on the electrical recording charts to assist in subsequent reading of the data. 

After the free-space calibration is completed, a static calibration is 
made on each sample material. The furnace is placed in the measurement 
position and each microwave sample in turn is rotated into the microwave beam. 

The variable precision phase shifter is adjusted to determine the room- temperature 
insertion phase for each sample. The transmission amplitude is determined by 
moving the sample axially in the beam through a distance of plus or minus one- 
quarter wavelength and finding the average transmission signal value. 

The natural-gas heating torches are lighted to begin the temperature 
cycle. The output of the control thermocouple is plotted on a recording chart, 
on which the desired temperature versus time profile has previously been drawn. 
Gas, air and oxygen flow rates are manually adjusted to follow as closely as 
possible the specified temperature profile curve. 

The following data are recorded versus time: optical temperatures, thermo- 

couple temperatures, change in sample insertion phase, reflection coefficient 
phase, sample insertion loss, incident power, and reflected power. During the 
run the sample is continuously moved axially in the microwave beam (± quarter- 
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wavelength longitudinal displacement) so that the average insertion loss can 
be determined. The electrical data are recorded on a six-channel visicorder, 
and the thermal data are recorded on separate recorders. An event mark signal 
is placed on all recordings once in each carrier revolution so that the 
response of each sample can be identified. A plot of the target temperature 
versus time profile is shown in Figure 7. 
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Figure 7. Plot of Ta 




DATA REDUCTION 


The electrical and the thermal analysis programs are written in American 
Standard Fortran for operation on a UNIVAC 1108 digital computer. This sec- 
tion serves to describe the theories and assumptions involved in the two 
computer programs. 

In order to make the required determinations of dielectric constant and 
loss tangents as functions of temperature, the electrical analysis program 
requires temperature profiles and sample thickness as a function of time for 
each material. The temperature profiles and sample thicknesses are obtained 
from a thermal analysis program. The temperature profiles and sample thick- 
nesses along with experimentally measured values of the power transmission 
coefficient and insertion phase delay are used in the electrical analysis 
computer program to aid in the determination of the piecewise approximations 
to the relative dielectric constant and loss tangent versus temperature curve 
for each dielectric material under consideration. This electrical analysis 
program is essentially the same as the one developed on two earlier Air Force 
programs (Ref. 1,2). Thermodynamic and transport physical property data 
with corresponding analytical equations, for all the sample materials are 
presented in this section. A listing of the thermal analysis program is 
given in the Appendix; the electrical programs have been published in Reference 
2 . 

Processing of Thermal Data 

The thermal analysis program is required to generate temperature versus 
distance profiles through the sample materials at various times during the 
experimental run. While the run is in progress, front surface, rear surface, 
and internal temperatures are recorded. The two surface temperatures and 
thermal transport properties of the materials are used in the thermal analysis 
computer program to calculate the required temperature profiles; internal 
measured temperatures are used to verify that the correct transport properties 
have been employed and that the program is functioning properly. 

Theory . In the measurements performed on this contract, the sample 
materials never exceeded temperatures on the order of 1530° K (2300° F) . 

Since they were considered candidates for the Space Shuttle thermal protection 
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system, it was reasonable to expect that no permanent densification or sub- 
stantial thermal expansion would occur. Thus, the thermal computer program 
was developed with the following assumptions; 

(1) Heat flow was one-dimensional in the axial direction, and 

(2) Density was constant and heat capacity and thermal conduc- 
tivity were functions of temperature. 

An energy balance over a thin slab of material whose flat faces are perpen- 
dicular to the direction of heat flow gives: 


pVc AT = [k(T - - T ) - k (T - T . ) 1 — 
K p n ' n+1 n n n_1 J Ax 


(5) 


where 

P 

c 

P 

k 

AT 

n 

At 

A 

V 

Ax 

n 


density of the material 
heat capacity 
thermal conductivity 

temperature in the slab, which we call the nth layer 
increment of time over which the heat balance is taken 
area of the slab perpendicular to the direction of heat flow 
volume of the slab 
thickness of the slab, and 

a distance index identifying successive slabs. 


Simplifying and noting that V = AAx, one obtains 


P c . 


AT - 
n 


kAT 

(Ax) 2 


(T + T - - 2T ) 
' n+1 n-1 n 


( 6 ) 


In this equation, AT^ represents the change in temperature of the nth slab, 
during the time interval At. Note that AT^ is uniquely expressed in terms of 
the thermal transport properties (p, c^ and k) , a time interval (At), the slab 
thickness (Ax) and three temperatures existing at the beginning of the time 
interval . 
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Various equations of this type could be written to describe the temper- 
ature in a slab, but the one shown was selected because it can be conveniently 
evaluated using a digital computer. The sample is divided into a number of 
imaginary slabs, usually 20. The transport properties of each slab are 
calculated based on its temperature. Then the temperature of each slab is 
calculated at the new value of time. The program continues to make time steps, 
calculating a new temperature versus distance profile after each time step. 

The front and rear surface temperatures are specified by input data. 

O'Brien, Hyman and Kaplan (Ref. 11) discuss the stability and convergence 
of numerical solutions of partial differential equations, and one of their 
examples is the one-dimensional heat conduction equation with a diffusivity 
(k/ PCp) of unity. They concluded that 10 to 20 layers are sufficient to 
achieve convergence (approach to the true solution), and also observed certain 
requirements for maintaining stability (damping of errors rather than growth 
of errors). The stability requirements impose limits on the relative sizes 
of time increments and slab thicknesses. 

Thermal Transport Properties of Samples . As stated earlier, density 
was assumed constant for all samples in this study; the fact that the 
temperature cycle was chosen to permit multiple cycling of the materials 
implied that no permanent change in density would occur. Densities used in 
the thermal analysis program are given in Table IV. 

i 

TABLE IV 

DENSITIES OF SAMPLE MATERIALS 


Material 

3 

Density (gm/cm ) 

SCFS 

1.902 

AS-3DX 

1.673 

Markite 3DQ 

1.913 

HD-0092 

1.976 

Mullite HCF 

0.279 

LI-1500 

0.237 
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The heat capacity of solids is a function of temperature and pressure 
only, provided the material does not experience a phase or structural change. 
Since all measurements on this program were made at a pressure of one atmosphere, 
the heat capacities were functions of temperature only. This dependence was 
represented by an equation of the form 

o 

c =a+bT+c/T 
P 

where 

Cp was expressed in units Btu/lb-m °R 
T was expressed in units °R, and 

a, b and c were constants derived from experimental data. 

The form of this equation results from considerations of the statistical 
thermodynamic description of solids; this subject is discussed in many standard 
thermodynamic texts such as Lewis and Randall (Ref. 12) . Experimental heat 
capacity data were furnished by NASA for LI-1500 and Mullite HCF. These data 
had been previously collected by Georgia Tech for the other materials (Ref. 2). 
Values of the heat capacity equation constants are given in Table V. 


TABLE V 

CONSTANTS FOR THE HEAT CAPACITY EQUATION* 


Material 

Constants 

a 

b 

c 

SCFS 

0,245 

0.142 x lO -4 

-0.250 x 10 5 

AS-3DX 

0.245 

0,142 x 10" 4 

-0.250 x 10 5 

Marki te 3DQ 

0.245 

0.142 x 10" 4 

-0.250 x 10 5 

HD- 009 2 

0.339 

0.514 x lO' 4 

-0.516 x 10 5 

Mullite HCF 

0.250 

0 

0 

LI-1500 

0.300 

0.810 x 10" 5 

| 

-0.480 x 10 5 


Heat capacity and temperature units are given in the text. 
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The thermal conductivity of amorphous or glassy materials has been 
reported by Jakob to be proportional to some power of the temperature (Ref. 13). 
Jakob also showed that for non -metallic crystalline substances the conductivity 
is inversely proportional to temperature. Materials investigated on this pro- 
gram included a combination of these types, so that a conductivity equation 
having this form was selected 

k = a + b/T + cT + dT 3 

where 

k was expressed in units Btu/hr ft °R 

T was expressed in units °R, and 

a, b, c and d were constants derived from experimental data. 

Experimental data were furnished by NASA for LI-1500 and Mullite HCF; they had 
previously been collected by Georgia Tech for other materials. Manual adjust- 
ment of conductivity was required in some cases to obtain satsifactory agree- 
ment among internal thermocouple temperatures and temperatures calculated by 
the thermal analysis computer program. The conductivity was chosen as the 
property to be adjusted in these cases, because it is generally less well 
established than density or heat capacity. For example, thermal conductivity 
is a function of sample density, and density often varies among individual 
specimens of the same material. Values of the conductivity equation constants 
are given in Table VI. 

TABLE VI 

CONSTANTS FOR THE THERMAL CONDUCTIVITY EQUATION’' 


Material 

Constants 

a 

b 

C 

d 

SCFS 

0.477 

0 

-0.183 x 10“ 4 

0.187 x 10 -10 

AS-3DX 

0.310 

0 

-0.119 x 10" 4 

0.122 x 10“ 10 

Markite 3DQ 

0.477 

0 

-0.183 x 10" 4 

0.187 x 10' 10 

HD- 0092 

9.060 

0.151 x 10 4 

0 

0 

Mullite HCF 

0.051 

0 

-0.167 x lO -4 

0.108 x 10" 10 

LI-1500 

0.078 

-0.207 x 10 2 

-0.181 x 10" 4 

0.106 x 10" 10 


«JU 

' Thermal conductivity and temperature units are given in the text. 
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Thermal Analysis Computer Program . The computer program receives input 
data in the form of front and rear surface temperatures at discrete times 
during the experimental run. Also, the appropriate thermal property constants 
are stored for use when needed. The sample is divided into a series of imag- 
inary layers (usually 20) and initial values of variables are established. 

The first time step is made (typically 0.5 second) and the two surface temper- 
atures at the new instant in time are calculated by linear interpolation of 
input times and temperatures. Then the temperature of each internal layer is 
calculated using equation (2) and the first profile is complete, A decision on 
whether to print the profile is made, another time step is taken, new surface 
temperatures are calculated, and the next set of internal temperatures is cal- 
culated; this completes the second profile. This process continues until the 
front surface begins final cooling (about 20 minutes after the start of the 
run), at which time the temperature profiles are not usable for reduction of 
electrical data. A flowchart of the thermal program is shown in Figure 8. 

Processing of Electrical Data 

Background . A homogeneous sample of high- temperature dielectric material 
is heated on one side so that the temperature distribution in the sample is 
given by T(x,t), where x is the distance from the "cold" side and t is time. 

The heated sample is positioned in a free-space , focused-beam waveguide bridge. 

A continuous record in time is made of the temperatures at discrete stations 
in the sample as well as of the measured values of power transmission coefficient 
and insertion phase delay of the sample. The duration of the measurements is 
on the order of 2400 seconds; i.e., the time required for the heated side of 
the sample to reach maximum temperature and cool to approximately 366° K (200° F) . 
It is desired to find the temperature dependence of the relative dielectric 
constant and loss tangent of the material from these measurements. 

A digital computer program has been prepared to aid in determining the 
piecewise linear approximations to the relative dielectric constant versus 
temperature curve and the log^Q (loss tangent) versus temperature curve for 
the material test sample. The program is written in Fortran IV for running on 
the Univac 1108 computer at Georgia Tech. 

Theory. An optimization technique is used to determine the piecewise 
linear approximations to the vs. T and log^Q (tan 6) vs. T curves which 
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Figure 8. Flow Chart for Thermal Analysis Computer Program 













minimizes the root mean square error between the measured and computed values 
of power transmission coefficient and insertion phase delay at discrete times, 
t : •••? in the experimental run. The vs. T curve is 

specified by the breakpoints of the linear segments comprising the curve. 

The temperature breakpoints (abscissas) are fixed by physical considerations; 
the ordinates of the breakpoints are allowed to vary. Thus, the vs. T 
curve is parameterized using a relatively small number of parameters jd: 

{ p l>P 2 > •••> P n 3 • For each sample, we choose n = 5. The log^Q (tan 6) vs. 

T curve is parameterized in a similar manner by the five parameters, £ : 

{ r l> r 2 > r 3 . For each sample, the same temperature breakpoints are 

chosen for the two curves. These parameters, £, r, are adjusted to bring about 
the optimization. 

The temperature at each of stations in the sample at each time, t^., 
is computed using the thermal program and the measured temperature data. The 
"cold" side of the sample is station number 1, the heated surface is station 
number The thickness of the sample at each time, t^, is also determined 


in the thermal program, 
program. 


These results serve as input data to the present 


For computation of insertion phase delay, , and power transmission 

c tC 

coefficient, c A^, at each time, t^., the sample is divided into 760 equal- 
thickness layers. The temperature at the midpoint of each layer is determined 
by orthogonal polynomial interpolation of the temperatures at the stations 
in the sample. (Subroutines ORTHLS and FITY.) The relative dielectric constant 
and loss tangent of each layer are then determined from the current vs. T 
and log^Q (tan 6) vs. T curves. The insertion phase delay and power trans- 
mission of the sample at time, t^, are then computed using the plane dielectric 
multilayers subroutine, MULTLY . 

From the above it is seen that the computed values of transmission co- 
efficient and insertion phase delay are functions of the parameters £ and _r; 
i . e . , 

A = c A k (E,r) (7) 


c k 


A 


(£.l) 


( 8 ) 


32 



Experience has shown that the insertion phase delay is a strong function of 
£ and has only a relatively weak dependence on r; also, the transmission 
coefficient depends more strongly on £ than on £. The optimization procedure 
is thus simplified by representing the computed values of insertion phase delay, 
e®, as functions of only £, and the as functions of only £. Since the vs. 

T curve has the dominant effect in these computations, optimization with respect 
to this curve is carried out before optimization with respect to the log^Q (tan 5) 
vs. T curve. 

The objectives of the computations are to determine the parameters £ and 
r such that 



number of discrete times, t— 

measured power transmission coefficient at t T . 

i\ 

measured insertion phase delay at t^ 
specified error tolerance in transmission coefficient 
specified error tolerance in insertion phase delay 

The error tolerances are based on measurement error. For the present version 
of the program, § A = 0.004 and 6x = 1.0 degree. 

A A A 

The parameters £ and £ which meet the requirements of (9) and (10) are 
determined in the following way. Initial estimates of the e r vs. T and log^Q 
(tan 6) vs. T curves are specified in terms of £ and £, and c A and are 
computed for these sets of parameters. The corrections Ap : {Ap^,Ap 2 , Ap^} 

such that 
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P 2 + Ap 2 ; 

» 

« 
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P n 


p + Ap 
n r n 


A _ 1 _ A 

= £ + Ael 


(ID 


are found using the following scheme. At a given time, t.., the differential 
of c \(£) is ’ by definition, 


d *. = 


a 5 K 
A c k 


A Pl + . . . + Sp 


8 }. 

c k 


n 


c k ' Sp 1 

The actual change in c $^(£.) for increments Ajj is 

A c $ k = c\ ( E- + A 2- } * c § k (£) 

For small A]>, we can equate (12) and (13) to yield 


Ap 


n 


( 12 ) 


(13) 


A § n - 


9 c\ 


S c\ 


k s ?1 Ap i + • • • + BPn 




(14) 


Let c \ ( £ + ^> = A 80 that 


Al k = m $ k ■ c 


(15) 


Doing this for each time t R results in the following matrix formulation: 


A$j| 


A . 

• • • # « 1 


" Ap l" 



9P]_ 

Sp n 



a$ 2 


• 

0 


Ap 2 

• 

• 

SB 

• 

* 

• 



» 

a§ LL 


0 

© 

34 LL 


Ap n 

_ 


dp-^ 

Bp n _ 




( 16 ) 
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ae = jf 1 M (19) 

where J ^ is the inverse of the approximate Jacobian matrix J. If LL > n, 
the system of equations is over-determined and the generalized unique inverse 
of the non-square matrix, J , must be used rather than the conventional inverse. 
Subroutine GID accomplished this. 

A 

Having obtained and, hence, £, the left side of equation (10) is com- 

A 

puted and compared to the specified error tolerance, 5^. In general, p) 

will not be equal to for all k = 1, 2, . LL. Hence, the process would 

A 

be repeated by finding a new set of corrections to £ found above. However, 
it is considered advantageous at this point to next apply the optimization 
procedure to the log^Q ( tan $) vs * T curve utilizing the r parameters and a 
duplication of the procedure outlined above for the £ parameters. In the actual 
computer program this is done; i.e. , the £ parameters are optimized first, then 
the £, then the £, etc., until equations (9) and (10) are satisfied. 

The incremental corrections, Ap , given by equation (19) would ideally 
be of the order of 5p since equation (14) and succeeding equations are valid 
only for small Ap . In addition, it is possible that equation (19) would yield 
values of Ap^ which would alter the vs. T curve to the point of being out- 
side physical bounds. Thus, in implementing this scheme on the computer, tests 
are provided to insure that the A£ are within physical bounds. In addition, 
a ’’fine adjustment” may be applied to the _A£ which consists of multiplying the 
Ap by a scalar constant, c, and adjusting this constant until 


A 

£ = £ + CAE. 


( 20 ) 
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yields a minimum RMS error in the insertion phase delay. Note that as c 
approaches zero, the corrected vs. T curve approaches the original curve. 

It was found during the course of the electrical data processing that no fine 
adjustment was required, resulting in a considerable savings in computational 
time. 

Computer Program . The computer program consists of a MAIN program and 
the following subroutines: FORWAR, MULTLY, GID, ORTHLS, FITY. The functions 

of these routines and the variable explanations are adequately described by 
the accompanying flow diagrams (Figures 9 and 10) and comment cards in the 
programs. Only a brief description of each of these programs is given below. 

MAIN is the calling program and carries out the optimization procedure 
described previously. FORWAR is called by MAIN to compute the power trans- 
mission coefficients, c A, and insertion phase delays, at the times, _t, 

for given sets of parameters £ and r. FORWAR accomplishes this by dividing 
the sample into NL = 760 layers, computing the temperature, relative dielectric 
constant, and loss tangent of each layer using ORTHLS and FITY, and computing 
c> A and of the resulting plane multilayer dielectric using MULTLY. MAIN 
builds the approximate Jacobian matrix, J, and utilizes GID to find the trans - 
pose of the generalized unique inverse of J. MAIN transposes LJ J to give 
J \ solves equation (19), and continues the optimization procedure until 
equations (9) and (10) are satisfied, or until other externally set conditions 
are met. 

The electrical program has been tested in the following way. A "true” 

e r vs. T curve and a "true" log^Q (tan 8) vs. T curve were assumed. These 

curves are piecewise linear curves specified by their abscissa and ordinate 

breakpoints. Typical input temperature and thickness data obtained from 

previous work were used as input test data. Using the "true" curves, the 

values of ^A were computed. These computed values were then defined as 

the "measured" data I, A. These values of $, A were then used as the 

m- 5 m— m- 5 m- 

measured input data. Rather gross estimates of the e and tan 6 curves were 
then assumed. These estimates consisted of linear continuations of the known 
low- temperature 255 to 1644° K (0 to 2500° F) , data into the high- temperature 
regions. The abscissa breakpoints on these estimated curves were the same 
as for the "true" curves. Using the same input temperature and thickness 
data as before, the program was allowed to run. The results show that only 
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Figure 9. Flow Chart for Electrical Analysis Computer Program. 
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Figure 10. Flow Chart for Subroutine Forwar 









one correction is required on each curve to satisfy equations (9) and (10) 
w hen 6^ = 1,0 degree and 6^ = 0.004, Setting these tolerances to zero and 
allowing the program to make one more correction to each curve resulted in a 
corrected curve essentially equal to the "true" curve. 

Processed Data . The measured data are surface and internal sample 
temperatures, insertion phase shift and transmission loss, and incident 
and reflected power. These are measured as functions of time for each tempera- 
ture cycle. The thermal data are then processed and the temperature profiles 
are obtained for each sample as a function of time and cycle number. 

The electrical data for each sample are plotted for each temperature 
cycle. Any deviations in electrical data from run to run are noted. The only 
significant changes occurred with hot-pressed boron nitride and the AS-3DX 
composite materials. Boron nitride lost material from the heated surface 
which amounted to approximately 1.6 electrical degrees in phase change per 
temperature cycle. The thickness of the boron nitride changed from 1 0 59 cm 
(0.627 inches) to 1.52 cm (0.600 inches) during the 10 temperature cycles or 
0.0068 cm (0.0027 inches) per cycle. The AS-3DX composite material insertion 
phase delay at room temperature was affected mostly by the initial temperature 
runs. The insertion phase delay decreased by 2 degrees after two runs were 
completed. The room temperature insertion phase delay then remained constant 
throughout the remainder of the test runs. 

The other materials (LI-1500, Mullite, Markite, and slip-cast fused 
silica) reacted in like manners. During the exposure to the heat, each material 
had phase shift and transmission loss responses that would be repeated from 
run to run. 

Typical data are plotted in Figures 11 through 16. These data indicate 
the transmission loss and insertion phase shift averaged over a number of runs. 
The run time indicated on each diagram can be correlated to temperature by the 
use of Figures 17 through 22, in which are shown the time- temperature history of 
the front surfaces of the samples. 

The remaining electrical data are presented as dielectric constant and 
loss tangent values as functions of sample temperature. As explained in the 
section on Processing of Electrical Data, the dielectric constant and loss 
tangent are obtained in the following manner. From the measurement run, sample 
temperature profiles are obtained at eight time intervals. The sample model 
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Figure 11. Georgia Tech SCFS: Carrier 1 and 4 Transmission Properties. 
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Figure 15. McDonnell Mullite HCF: Carriers 2 and 3 Transmission 

Properties . 
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is divided into 760 layers for each time interval. Assumed dielectric constant 
and loss tangent values as functions of temperature for the sample are used 
and each layer is assigned specific values of and tan 6. The multi-layer 
analysis transmission computer program is then used as a tool to calculate the 
insertion phase shift and transmission coefficient for these conditions. These 
calculated data are then compared to the measured insertion phase shift and 
transmission loss. An error term is generated which is used to perturb the 
assumed dielectric constant and loss tangent versus temperature curves. 
Iterations are made until the calculated data approach the measured data. In 
this manner, the end results are curves of dielectric constant and loss tangent 
versus temperature for each measured sample. The short-circuit waveguide data 
obtained from NASA Langley are used as a baseline in the determination of the 
initial temperature functions. The computed values of dielectric constant and 
loss tangent are presented in the section entitled Results and Discussion. 
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RESULTS AND DISCUSSION 


Four carrier disks were run for collection of experimental data. Each 
carrier was exposed to the specified reentry temperature profile ten times, 
using the temperature of the slip cast fused silica surface for control pur- 
poses. (When a carrier contained no SCFS sample, the carrier surface tempera- 
ture was used for control.) Carrier and sample assignments are given earlier 
in Table III. 

The purpose of this program was to determine whether a deterioration of 
dielectric properties occurred in the candidate window materials during re- 
peated recycling through a simulated reentry temperature versus time profile. 
During a typical experimental run, data were recorded for about 18 minutes; 
there were a total of 40 experimental runs, each involving three sample mate- 
rials. In view of the magnitude of the data reduction task, it was decided 
to begin by processing only the first and tenth runs of each carrier. If no 
deterioration of dielectric properties occurred between the first and tenth 
runs, then no further information would be gained by processing the inter- 
mediate runs. If deterioration was observed, then intermediate runs could be 
processed as required. Runs 1, 5 and 10 of the first carrier were evaluated, 
and Runs 1 and 10 of subsequent carriers were completely processed. 

Thermal Data 

Figures 17 through 22 show typical temperature versus time profiles for 
each of the sample materials. The three silica based materials (SCFS, AS-3DX 
and Markite 3DQ) all showed the same temperature response behavior to the 
heating environment. We were unable to reach the rate of rise desired at the 
beginning of the mins, but successfully matched the desired maximum tempera- 
tures. We believe that the temperature exposure was a reasonable simulation 
of the desired reentry profile, however. Note also that the difference between 
front and back surface temperatures was on the order of 200° K (400° F) through- 
out the runs for these materials. 

The boron nitride (HD-0092) samples experienced a somewhat slower rate of 
front surface temperature rise than did the silicas, but they also reached the 
desired maximum temperature. The difference between front and rear surface 
temperatures was on the order of 100° K (200° F) . The lower rate of rise and 
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Figure 17. Typical Temperature versus Time Plot for Slip-Cast 
Fused Silica. 
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Figure 18. Typical Temperature versus Time Plot for AS-3DX. 
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Figure 19. Typical Temperature versus Time Plot for Markite 3DQ. 
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Figure 21. Typical Temperature versus Time Plot for Mullite HCF 








the smaller temperature gradient can both be attributed to the very high ther- 
mal conductivity of boron nitride. 

The low density samples (Mullite HCF and LI-1500) both experienced high 
rates of front surface temperature rise, almost matching the target tempera- 
ture profile. They reached the desired maximum temperature and had steep 
temperature gradients, on the order of 400° K (700° F) between the front and 
rear surfaces. The steep gradients and high rate of temperature rise can be 
attributed to their low thermal conductivity. 

These results illustrate the influence of the thermal properties of 
materials on their temperature response to a reentry environment. Various 
materials responded differently even though heating conditions were reproduced 
as closely as possible; in fact different materials run on the same carrier 
disk showed different rates of temperature rise. There is no "representative 
dielectric window material" for thermal response calculations. Reentry temper- 
ature calculations must be run using appropriate thermal properties for each 
candidate material, to accurately predict the temperature which will be reached 
during a proposed reentry flight path. 

Of the six specimen materials, only two showed any physical evidence of 
change after ten temperature cycles. Boron nitride (HD-0092) had a powdery 
deposit on both faces after removal from the carrier disk. The deposit on 
the front surface showed streaks, and was an estimated 0.5 mm thick; the de- 
posit on the rear surface was smooth and somewhat thinner 0 Mullite HCF had 
lost some flakes of coating around the perimeter of the sample after removal 
from the carrier, and the remaining coating was loosely held to the porous sub- 
strate. Figure 23 shows microwave samples of each material after ten tempera- 
ture cycle So 

Electrical Data 

The results of the microwave transmission measurements are presented in 
this section as dielectric constant and loss tangent of each sample versus 
temperature. The dielectric constant and loss tangent of McDonnell Mullite 
HCF vary with temperature as indicated in Figure 24. The dielectric constant 
increases by approximately 7 percent from room temperature to 1480° K (2200° F) . 
The loss tangent increases from a room temperature value of 0.0003 to 0.068 at 
1480° K (2200° F) o This is a large change in loss tangent as compared to the 
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Figure 23. Sample Materials After Ten Temperature Cycles 
(from upper left: SCFS, AS-3DX, Markite 3DQ , 

Boron Nitride HD-0092, Mullite HCF, LI-1500) . 
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other materials considered. 

In Figure 25 are illustrated the electrical properties of Lockheed LI- 1500. 
The dielectric constant increased by approximately 2 percent from room tempera- 
ture to 1530° K (2300° F) . The loss tangent increased from 0.0005 at 298° K 
(75° F) to 0.0015 at 1530° K (2300° F) . Thus, the loss tangent increased very 
little over this temperature range. 

The electrical properties of General Electric Markite are presented in 
Figure 26. The dielectric constant increases almost linearly from a 298° K 
(75° F) value of 3.17 to a value of 3.26 at 1255° K (1800° F) . The loss tan- 
gent shows no increase with temperature to 1255° K (1800° F) . Above 1255° K 
(1800° F) , both the dielectric constant and loss tangent increase. The loss 
tangent shows a rapid increase. 

The Philco-Ford AS-3DX material shows little change in either dielectric 
constant or loss tangent over the measured temperature range. These data are 
presented in Figure 27. 

The electrical properties of the Georgia Tech slip-cast fused silica are 
illustrated in Figure 28. The dielectric constant increases only slightly up 
to 1366° K (2000° F) . Above 1366° K (2000° F) , the dielectric content in- 
creases at a faster rate with temperature. The loss tangent indicates a linear 
Increase from a value of 0.0005 at 298° K (75° F) to a value of 0.006 at 1644° 

K (2500° F). 

In Figure 29 are presented the electrical properties of hot-pressed boron 
nitride. As indicated, the dielectric constant remains almost constant 
(e = 3.93) to 810° K (1000° F) and from that point, the dielectric constant 
increases to a value of 4.05 at 1394° K (2050° F). The loss tangent remains 
very low from 298° K (75° F) to 588° K (600° F) . It increases at a linear 
rate to 1394° K (2050° F) . As previously stated, the boron nitride did ablate 
slightly during each measurement run. 

All of the materials except McDonnell Mullite exhibit excellent low loss 
properties as functions of temperature. The McDonnell Mullite had a much 
higher loss tangent than the other materials tested, but this does not prohibit 
its use as a window material. This is especially true if other factors make 
it a better candidate for a particular application. 
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CONCLUSIONS 


The following conclusions seem warranted from the results of this measure- 
ment program: 

(1) A system was developed which permitted measurement of the micro- 
wave properties of candidate window materials during simulation of the Space 
Shuttle reentry thermal profile. 

(2) None of the six materials tested showed significant deteriora- 
tion in dielectric properties with repeated cycling through the specified re- 
entry temperature profile, up to ten temperature cycles. AS-3DX showed small 
changes during the first two cycles. 

(3) Four of the sample materials showed no visible deterioration in 
physical properties after ten reentry temperature cycles. Mullite HCF suffered 
loosening and chipping of its emittance control coating which might be serious 
in the presence of aerodynamic drag. Boron nitride lost thickness and devel- 
oped powdery surface deposits, probably the result of oxidation; however, this 
does not appear to be a severe problem after ten temperature cycles. 

(4) The dielectric constants and loss tangents of the sample mate- 
rials as functions of temperature are within ranges that would not prohibit 
their use for Space Shuttle antenna windows. Mullite HCF had a higher loss 
tangent than the other candidates. 

(5) A smaller program would be advisable to verify the acceptability 
of the final dielectric window design, once the material and reentry tempera- 
ture profile have been firmly established. 
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APPENDIX 


Thermal Analysis Computer Program 





oom 


00106 


00106 


00106 


ll FORMAT (lHl) 

- la F0 rmAT(/ > 6X>11 HM ATE R IAl = .9A6 ./ *-9 X> 12 H T H ICKNESS = , F5 t 3>5 H INCH) 
13 FORMAT (/>9Xr 26HHEAT CAPACITY C0EFFlCIENTS./>l2X*4HA = »El0.3»/» 


20EFFICIENtS,/»12X»4HA : 


49HDELTIM s ,F6,3»4H SEC»/»9 Xp23HNUM0ER OF TIME STEPS = »I5» 


00107 


UUAiO 

QOlll 


"UUXXX 

00112 


yu* xc 

00113 


00115 


AW 

00117 


00120 


00130 


111 F0RMAT(l3X»F7.1*14XfF6.1*17X»F6,l) 

. _ — m t , CO ■» . ii - tj — g »gf> / \ - 


1 ft FORMAT (/F9Xf7HSTATI0Nf6X»13HDlSTANcE ( IN) • 6Xf 12HTEMP (DEG F>#6Xf 

. <u«. . on g-U*y* I gfl I — 


17 FORMAT tl2X > I2 > llX»F5.3>l‘fX>F5.0»15x>E9.3>14X»F6.‘U 


2l FORMAT (9A6»F10» 3) 

~ 2g FORMAT (3E ~ 10. 31 ~~ 

23 FORMAT (4E.10.3) 

- 24 FORMAT (2F7»3 ) 

25 FOrmAT(3F10,0) 


READ INPUT DATA , 

READ (5.22) CArC0.CC 



00225 

00227 

00232 


DI ST ( I ) = ( DELX*FUOAt ( I ) -DELX ) *12. 0 
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THERMAL ANALYSIS PROGRAM (Continued) 


—00233 

00235 

58* 

59* 


"3sr 

TtrASm) = 540.0 — — — 

TIM = DELTIM 

00236 

00237 

6tr* 

61* 



_jg ■ ■ ■ ■ 

J3 : 0 

0024Q 

00241 

62* 

63* 



"t*0 

TIM1=0.0 

00242 

00243 

64* ~~ 
65* 



"TIM2 = 0tO 

TS2=80.0 

00244 

00244 

66* 

67* 

c 


Tli£«»uo*Q 
START time loop 

00245 

00245 

66* 

69* 

c 


“Dv 59 J=1 9 “ ■ • • ■ — - ~ ■ 

calculate front and back surface temps, when tim has passed the 

0024b 

00250 

y.Q-jU 

71* 

c 

f» 


"END OF PRESENT INTERPOLATION LINE? CALCULATE NEW" LINE — 

IF { Tl M-TI M2 ) 35 » 35 * 34 

qq25u 

00253 

72* 

73* 


34 

■c«l^.CUt:RTE "Ew INTtKrOLATION LINE — 

TSi=TS2 

“00254 

74* 



"Toi=:TB2 

00255 — 
00255 

75* 

76* 

c 


READ a new data point 

0025b 

00257 

77 * 

78* 



“L = L+1 

TlM2=TIML(L) 

— 00260 
00261 

79* 

80 * 



"TS 2 -TSI L > ” “ ' ' ' 

TB 2 =T 8 (L) 

00262 

00263 

81* 

82* 



P5lO“E5(T52-TS1)/*TIM2-TX'M1) _ 
BSlOPEs ( TB 2-TB1 ) / ( TIM2-TI Ml ) 

00264 

00265 

26* 

84* 



FInTCPbTS1"FSLOPE*TIM1+480. 0 
8IntCP=TB1-3SLOPE*TIM1+460.0 

00266 

85* 

86 * 

c 

35 

'FIND SURFACE TEMP FROM LINE 
T(21)=FSL0PE*TIM+FINTCP 

— 00267 — 
00267 

_ — $7* 

88 * 

c 


T 1 1) =BSLOPE#TIM*BImCP 
CALCULATE INTERIOR STATION TEMPS 

— 0027TJ 

00273 

89* 

90* 


•Jf JU 

K=KA + < KB/TLAST ( I ) ) +KC*TLAST ( I ) +K0* ( TL AST ( I ) **3 ) 

00275 

91* 

92* 



ALPHA <I1=K/<CP*RH0) 

0027b 

0027? 

98* 

94* 


*0 

RftLr A l J ) =R*AtrHA IT J ' “ " 

T ( I ) =TLAST ( I ) +ALPHA (I ) *R* ( TLAST .( I +1 ) +TL AST ( 1-1 ) -2*TL*ST ( I ) ) 

00277 — 

00301 

95* 

96* 

- c- 


“DECIDE WHETHER TO PRINT* WHEN J 2 = M 2 , PRINT 
IF(J2-M2), 42,43,43 

00304 — 
00305 

97* 

98* 



J 2 sj 2 +l 

SO TO 48 

00303 

00306 

— 99* — 
100* 

c 

*3 

it 

UC^ xut. — writ * nfcK — I U "l UtlPI “not f 1 f 1 U“'In r RUE- 

XF( J3«i ) 4^,45,45 

00311 

00313 

— tor* — 
102 * 


9tr 

WKI 'OF 11* 

J3=l 

UU314 

00315 

ft ft H ^ 

1 UO* 

104* 

1 ACa. 


*5 

11 ^ 

“tfV IV 

J3 = 0 

i.lOTtr/7 .TCi «i»T*A 

U U Ol & 
00321 

X U D* 
106* 



1 & * 0 9 XZt * | l! ¥ J — 

WRITE (6, 16) 

' 00323 
00326 

1 U / ^ 
108* 



~£rv — 4 *t J-S-l f &X ! 1 

TLaST(I)=T(I) 

00327 

00330 

X U 7* 
110 * 

/•* 

4? 

T^IJ—T ( 1 J j **+qU • O 

WRITE (6,17)1,01 ST (I) ,T(I) *ALPHA(I) ,RALfA(I) 

nrrt r» 1 lM^.U etilAl QOftPTI re APlH UlOrfF AklD ffftDUAT 0*7 f*Ap»i-wC UETPJT 

00*530 

00340 

ill* 

112 * 

L 


J 2 =l 

00041 

00341 

li3* 

114* 

c 

/> - 

sy ry so 

SET up for next time loop 

00342 

00345 

11 a* 

116* 


**8 

49 

"t)v 1+7 X — 1 9 Cl 

TLasT(I)=T(I) 

— 00347 — 
00351 

117* 

118* 


53 0 

Tiiis i JrMtytti rtFl 1 

NR UN = NRUN+1 

00351 

00352 

tt9* 

120 * 

c 


"WHEN NRUN=NRUN5» TOU ARE" FINISHED 
IF(NRUN-NRUNS) 30*60,60 

00357 

— 12 t* — 
122 * 


— &xr 

WRITE(6,11) 

STOP 

0035 1 
00360 

123* 

124* 

c 


"DATA CARDS FOLLOW XQT CARD 

END 


END OF COMPILATION*. 


NO DIAGNOSTICS. 
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